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СПОСІБ ВІДНОВЛЕННЯ ПОТЕНЦІАЛУ 
ЗА ЗНАЧЕННЯМИ МОДУЛЯ ЙОГО ГРАДІЄНТА 


Однією з важливих задач теорії потенціалу, що виникають при інтерпретації даних 
гравіметрії та в теорії фігури Землі, є гранична задача відновлення потенціалу сили тяжіння 
за значеннями модуля його градієнта, якими, по суті, є виміряні значення аномалій сили тя- 
жіння. Теоретичне підгрунтя цієї задачі закладено на основі способу аналітичного продов- 
ження модуля градієнта потенціалу сили тяжіння (МГІПСТ,) для (лінійного) диференціально- 
го рівняння сили тяжіння, виведеного А.В. Чорним |1|. Ця гранична задача має два альтерна- 
тивні шляхи розв'язання: ітераційний пошук границі послідовності розв'язків лінійних задач 
Неймана для рівняння Лапласа, реалізований А.І. Якимчиком |2|, та постановка так званої 
нелінійної задачі Алексідзе для того ж рівняння на довільних банахових просторах, що є ме- 
тою даної замітки. 

Відомо, що потенціал сили тяжіння У(х) є гармонічною функцією в області С", не 
зайнятій тяжіючими масами Землі. Однак скористатись цією його властивістю для відновле- 
ння значень сили тяжіння «(х) в певних точках х необмеженої області С" шляхом розв'яза- 


ння зовнішньої граничної задачі Неймана для рівняння Лапласа (це стосується так само і змі- 


з . ду (х) 
шаної задачі) з граничною умовою ант 2(х), х є ДУМИ, заважає неспівпадіння земного ре- 
п 


льєфу СС - поверхні вимірювань аномалій «(х) з будь-якою еквіпотенціальною поверхнею 
ОМ: У(2)- Сх, тобто, наближеність граничної умови при точному диференціальному рів- 
нянні спричиняє неточний розв'язок задачі Неймана. Так само непридатний для редукції зна- 
чень «(х) з поверхні СС в область С" і розв'язок зовнішньої задачі Діріхле для рівняння 
Лапласа (заодно 1 складнішої задачі Коші для того ж рівняння) через негармонічність функції 
а (х), х є С", відзначену М.А. Алексідзе |31, тобто маємо наближене диференціальне рівнян- 
ня і точні граничні дані. 

Запропонована К. Юнгом |4| підміна приросту МГПСТ в локальній області вертика- 
льною складовою сили тяжіння дозволила використати в тлумаченні аномалій потенціалу 
апарат теорії гармонічних функцій. Така трактовка значень аномалій де(х) як значень гар- 
монічної функції (в задачах гравірозвідки), її нормальної похідної, або їх лінійної комбінації 
(в задачах геодезичної гравіметрії), що задані на частині чи всій границі області її визначен- 
ня, одначе, при розв'язанні граничних задач справедлива лише для областей малої міри 
(19х1? з глибиною до 10 км, сфера гравірозвідки, див. (21). Перенесення такої підміни на по- 
ширені регіональні гравіметричні побудови, тим паче з метою вивчення розподілу глибин- 
них неоднорідностей в корі і верхній мантії, неправомірне: відхилення від істинного розв'яз- 
ку прямо пропорційне геометричним параметрам області, що вивчається. 

В роботі |3| вперше з часів Юнга М. Алексідзе нагадав, що аномалії Де (х) є значення- 
ми МГИПСТ, а не однією із складових цього градієнта та вказав наближений спосіб відшукан- 
ня гармонічного наближення до реального розподілу де(х) в локальній області з точністю 
до сталої та геометрії самої області. Оскільки такі сталі залежать від форми й розмірів облас- 
ті малої міри (істотно різняться для рівнин 1 гірської місцевості), то в разі відновлення Де (х) 
в глобальній області з її значень в локальних підобластях відсутні будь-які критерії задовіль- 
ного ,, склеювання" гармонічних наближень і відповідні побудови стають необгрунтованими. 

Відтак, на основі даних про де(х) неможливо сформулювати точні граничні умови 


для жодної з класичних граничних задач математичної фізики (в т.ч. і для задачі Молоденсь- 
кого). Для характеристики розподілу да(х) в глобальній області слід повніше враховувати 
диференціальні властивості сили тяжіння. Узагальнюючи вищевказані проблеми редукції 
значень Ля (х), фактично отримуємо таку альтернативу: 


1. віднайти диференціальне рівняння, якому б задовольняли значення МГПСТ, і вирішити 
для нього відповідну граничну задачу; 
2. вирішити таку граничну задачу для потенціалу сили тяжіння, в граничних умовах якої 


безпосередньо закладено значення 2(х) хєС". 


Перший варіант цієї альтернативи знайшов своє втілення у роботі |2|, а наразі на загал 
пропонуються міркування щодо другого варіанту вирішення вказаної проблеми. 


Нехай С" - обмежена область точок тривимірного евклідового простору, зайнята ма- 
сами Землі, С" - необмежене доповнення до неї, вільне від тяжіючих об'єктів, СС - границя 
між ними, тотожна фізичній поверхні Землі. Маси М (2) Сє С" всередині Землі густини 
ам (2) - с(бмє генерують зовні неї поле сили тяжіння з потенціалом 


Ох) хеС 


0, хеєт" 


сія) 
(а тає а) 
я їх -б | 
У 
де 7 - гравітаційна стала, О(х)- 0.5а7 (х, - х, ) - потенціал центробіжної сили, шФ - модуль 
вектора кутової швидкості Землі, при цьому потенціал У()є С І, (сг) належить до класу 
неперервних функцій разом з похідними до другого порядку включно, якщо густина тяжію- 
чих мас с(2) є С сі ) належить до подібного класу разом з першими похідними. 
Напруженість поля (значення модуля градієнта потенціалу) дорівнює 
3, СУ (х) дж, (х) 
- 22 ан еВ 
32 УМХ) -У, -(з)л()), 2) 
ее дх, дп 
дж. (х) і ре і 
де око щ сок(п, ху , Кк - 1,2,3 - напрямні косинуси одиничного вектора л(х) внутрішньої но- 
п 
рмалі до еквіпотенціальної поверхні ОМ, : М.) - Сх, яка проходить через точку х. 
В рамках цієї моделі МГПСТ задовольняє (нелінійному в загальному випадку) дифе- 
ренціальному рівнянню сили тяжіння (еліптичне другого порядку типу Гельмгольца): 


у б'ю) Р В - 4д|ягад с(х), хе С" З) 


2 
стодх 0, хеСс" 
де я 5) - геометрична характеристика (міра кривизни) поверхні СИ : Ух) Сх. Саме це 
рівняння лягло в основу вирішення першого варіанту альтернативи: знайти функцію «(х), 


х є С" , яка задовольняє всередині необмеженої замкненої області С" рівнянню (3) за умови 
Ши «(х) -» 0 за її значеннями «(х) С (2С) , хе дС вбудь-якій точці границі СС. 


хо 
Альтернативний підхід ілюструє така постановка задачі: знайти функцію МУ), хеС", 
яка задовольняє всередині необмеженої замкнутої області С" - С" СС рівнянню Лапласа 
ДУ (ж) - 0, х є С", якщо в будь-якій точці ляпуновської границі СС області і в нескінченно 
віддаленій точці вона задовольняє умови 
3 2 
У ото) - ах), хе дО, У(х)-» 0, якщо |» о, (4) 
ст дх, 
де х(х) - задана неперервна функція. Назвемо цю нелінійну граничну задачу задачею Алек- 
сідзе для рівняння Лапласа. 
Відшукуючи гармонічну в області С" функцію Уу(х) у вигляді потенціалу простого 





шару У(х) з - оіб) 45., хє С" невідомої густини с(х), хедС в результаті серії аналі- 
Ал 4, 

тичних перетворень отримуємо щодо шуканої густини наступне нелінійне інтегральне рів- 

няння сили тяжіння: 





ПОР с(х) со5(т, р) 
, (х)-н ри і ут с(2)45. - 


с(п) 2 
- ,хедс 5 
раю Й. Я й , со5(р, 415 -45, 2 (х) хє (09) 


Розв'язок цього інтегрального рівняння еквівалентний розв'язанню задачі Алексідзе з 
граничними даними на поверхні Ляпунова, оскільки за будь-якого вибору густини потенціал 


Сх) простого шару задовольняє в С" рівнянню Лапласа, а знайдене з останнього рівняння 





сили тяжіння значення густини забезпечує виконання граничної умови (4). Питання розв'яз- 
ності задачі Алексідзе цим зводиться до вияснення умов коректності рівняння (3). На основі 
властивостей й кро б рівняння (5) можна спростити до вигляду 


о | ге а со5(р, а)45, 45; -28"(х) хедС. (6) 


дсдеХ 1 





За цими двома рівняннями можливо знаходити не лише значення потенціалу Ух), але й зна- 


чення модуля його градієнта в будь-якій точці області С", які можна обчислити з виразу 


згода У кое м, |, 


167" Кеї РАРАЄ - 


який не потребує обчислення чотирикратного інтегралу, як у випадку рівняння (4). 
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